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Introduction 

Cauchy problems for weakly dissipative wave equations 

utt-Au + b{t)ut = 0, u{0,-)=ui, Btu{0,-)=U2, (0.1) 

where the coefficient b = b{t) is assumed to be positive and tends to zero as t tends to 
infinity, provide an important model problem for the study of asymptotic behaviours 
and the influence of lower order terms on them. We refer to |RW04j for an exposition 
of results in that case or to the papers of A. Matsumura |Mat77j . H. Uesaka |Ues8flj or 
F. Hirosawa and H. Nakazawa |HN03j for known energy decay estimates. 

In special cases solution representations are explicitly known and thus also the question 
of sharpness of estimates can be answered. For properties of free waves we refer to the 
book of R. Racke |Rac92j or the classical paper of R.S. Strichartz, |Str70j . On the other 
hand the case of scale invariant weak dissipation, b{t) ~ t~^ as t ^ oo, is closely related 
to Bessel's differential equation and considered by the author in |Wir04j . 

The purpose of this paper is to prove structural properties of representations of solutions 
in the case below scaling, i.e. we assume that the coefficient function decays faster than 
t~^. This case is closely related to free waves. In case of integrable coefficients it is known 
from the works of K. Mochizuki and coauthors |Moc77j . |MN96j or from the preprint 
|Wir02j that the solutions are asymptotically free and a scattering theory in the sense 
P.D. Lax and R.S. Philipps |LP73j can be obtained. 

We will employ the translation invariance of the Cauchy problem. This implies that a 
partial Fourier transform with respect to the spatial variables may be used to reduce 
the partial differential equation in u{t,x) to an ordinary differential equation for u{t,^) 
parameterised by the frequency parameter |^|, 

uu + b{t)ut + \(\\ = 0. 

Its solution can be represented in the form 

u{t,0 = Mt,0ui + <^2it,0u2 

in terms of the Cauchy data ui and U2 with suitable functions (Fourier multiplier) $1 
and $2- 

In general estimates for $ are complicated to obtain, so the natural starting point is 
to rewrite the second order equation as system for the energy {\^\u,Dtu)'^ or a related 
vector and to use a diagonalization technique to simplify the structure and to estimate 
its fundamental solution. 

Main results of this paper are the solution representation of Theorem II . 101 together with 
its consequences for the U'-L'^ decay. Theorem 12.51 Furthermore the sharpness of these 
results follows from a modified scattering theory given by Theorem 13.21 
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1 Representation of Solutions 

We start by giving examples. The technique used in this paper allows to consider the 
following coefficient functions. 

Example 1.1. Scale invariant dissipation, b{t) = j^ is contained for fx £ [0, 1). For 
the general case /i > we refer to |Wir04j . The exceptional value /x = 1 is critical 
for estimates of the solution itself and related to an exceptional behaviour of Bessel 
functions. 

Example 1.2. Let /i > and n > 1. Then we can consider 



m 



(l + i)log(e + t)---logl"J(eN+t) 



Here we denote by log'"^ the n times iterated logarithm, log^ '(r) = r and log''^ '{t) = 
log log' '{t), and by e'-^' similar the iterated exponential e'*^^ = 1 and e^ ~^^' = e^^ ' . 

Exam,ple 1.3. Monotonicity of the coefficient is not essential. We can also consider 

coefficient functions like 

^ 2 + cos(alog(e + t)) 

^' 4(e + t) 

with a e M. 



Assumptions. We make the following assumptions on the coefficient function h = b{t) 

(Al) positivity: b{t) > 0, 

(A2) symbol-like estimates: for all A; G N it holds 



at^'^'^ 



<Ck 



I \ l+A; 



,1 + i. 

(A3) it holds \mlsn]i^^^tb{t) < 1. 
The last assumption is necessary to exclude the above mentioned critical case b{t) = jr^. 

Basic ideas, zones. If we consider the vector valued function U = {\^\u,'Dtu) , the 
energy, we obtain the system 

In order to understand the properties of its solutions we relate the size of b{t) with the 
size of \^\. This leads to the introduction of different zones in the extended phase space. 
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Therefore, we will use the implicitly defined function 

{i+t^m = N (1.2) 

with suitable constant N and define the 

ZhypiN) = {it,0\t>t^}, Zdiss{N) = {{t,C)\0<t<t^}. (1.3) 

Inside Zdiss the positivity of the dissipation b{t) is essential and we will use a reformu- 
lation as integral equation to solve the system. 

In Zfiyp the non-diagonal entries are dominating and we use a diagonalization scheme 
to transform the system to a diagonally dominated one. In order to do this we employ 
the symbol estimate of Assumption (A2) and introduce symbol classes in the hyperbolic 
zone Zhyp- This is related to the approach in lEYOflj . 

Definition 1.1. The time- dependent Fourier multiplier a{t,£,) belongs to the hyperbolic 
symbol class SN{mi,m2} if it is supported in the hyperbolic zone and satisfies the symbol 
estimates 

/ 1 \ m.2+k 

for all {t, G Zhyp{N) and all ken, a e N". 

Thus using (|1.2() we can conclude the essential embedding rule 

SNi'mi — k,m2 + i} ^^ SN{mi,m2} \/ i > k > 0. (1-5) 

Definition II . II extents immediately to matrix- valued Fourier multiplier. The rules of the 
symbolic calculus are obvious and collected in the following proposition. 

Proposition 1.1. 1. SN{mi,m2} is a vector space. 

2. 5Ar{?ni, 7712} • SN{m[,m'2} ^^ S^lmi + m[, m2 + 1^,2} 

3. DjD|S'Ar{mi,m2} ^^ S'7v{mi - |a|,m2 + A;} 

4. SN{-l,2]^LfL] 



For later use we define the larger symbol classes 5"^' ^{fni, 7712} of symbols with restricted 
smoothness. For these symbols we assume the symbol estimate H1.4() for k = 0, 1, ... £1 
and \a\ < £2- These symbols with restricted smoothness can be used to deduce map- 
ping properties in L^ spaces^. We give one auxiliary result following directly from 
Marcinkiewicz multiplier theorem, |Ste7n| Chapter IV. 3, Theorem 3]. 

Proposition 1.2. Each a G 5"^ ^ {0,m} gives rise to an operator a(f , D) : L^ ^ L^ 

for all p £ (1,00) with norm estimate 

||a(t,D)||p^p< r^ 



'^We are speaking about Fourier multiplier only, so no essential difficulties can arise by this lack of 
smoothness. 
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Formulation in system form. We do not use the vector function U mentioned before. 
The main reason is that we prefer a representation, where we can extract also estimates 
for the solution itself. Furthermore for small frequencies we need information on the size 
of the solution in order to obtain sharp decay results. 

We consider the micro-energy 

U = {h{t,Ou,Dtuf (1.6) 

with 

N 
^(*' ^) = Y+t'^'^*'"'^^*' ^^ ^ \^\'l'hyp,N(.t, 0- (1-7) 

Here and thereafter we denote by (l)diss,Nit,0 the characteristic function of the dissipa- 
tive zone and by 4>hyp,N{t,0 the characteristic function of the hyperbolic zone. 

This micro-energy is related to the energy operator E(t,D), 

E(t, D) : ((D) ui,U2f ^ mu{t, ■),Btuit, ■)f, (1.8) 

also considered in |Wirn4j . For all t we estimate some kind of H^-norm., but uniformly 
for all t the estimate can only be seen as H^ estimate. 

The aim is to prove estimates and structural properties for the fundamental solution 
S{t,s,^) to the corresponding system DjC/ = A{t,S,)U. 

1.1 The dissipative zone 

The essential idea is to write the problem as an integral equation. In the dissipative 
zone the micro-energy H1.6() simplifies to 

/ N 

U = u,Dtu 

\l + t ' 



and thus we have to solve the system 



i N \ 

Bt£{t, s, = A{t, 0£{t, s,0={ (i^fP ,^+*, 1 £{t, s, 0, £{s, s, = I. (1.9) 

N 



(i±^ ^b{t)] 



We rewrite it as an integral equation. Let therefore A(t, s) be the fundamental solution 
of the diagonal part of this system, 



'1; 

i+t 



A{t,s)= '^' xHs) , (1.10) 

X{t) = exp s ^ /q 6(r)dr >, such that the application of Duhamel's principle yields 

£{t, s, = A(t, s)+ij A(t, T)RiT, OS{t, s, e)dr (1.11) 
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with R{t,^) = A{t,$^) — diag A{t,^). Assumption (Al) guarantees that the entries of 
A(t, s) are monotonous functions. 

We claim that the dominating entry of A(t, s) describes the behaviour of £{t, s, ^). If we 
consider A{s,t)£{t,s,^) we get the integral equation 

A(s, t)£{t, s,0=I + i f A{s,t)R{ t ,OA{t,sI (A(s, t)£{t, s, 0) dr (1.12) 

with kernel K{t,s,£,). We will prove well-posedness of this equation in L°° over the 
dissipative zone using Theorem lA.il It holds 



A!(T)(l+£)|t|2 



^(^^^^e)=l AV)(i+.) „„ ^"^'^^]. (1.13) 



Therefore we have to relate the auxiliary function A^(t) to t. This is done using Assump- 
tion (A3) and the following consequence of it. 

Proposition 1.3. Assumptions (Al), (A3) imply 

dr t 



/o A2(r) A2(t) 
and jTTf) is monotonous increasing for large t. 

Proof. Integration by parts yields 

* dr _ _^ /■* TbJT) 

On the one hand the right hand side is larger than tA~^(t) by Assumption (Al). On the 
other hand we conclude from tb{t) < c < 1 for t > to that 

dr < / -T^dr + c / -r—- <C + c ' 



/o A2(r) - 7o A2(r) ^ A2(r) " 7o A2(r) 

and the statement follows 

Jo xHr)-i-c[ ^ x^it))- x^ity 

For small t the statement can be concluded from A^(t) ~ 1. 
Monotonicity is a consequence of 

d t _ 1 - tb{t) 
dt A2(i) ~ A2(t) 

and i6(t) < 1 for t > 1. D 
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Especially under Condition (A3) we have A^(t) < (1 +t). 
Lemma 1.4. Assume (Al) and (A3). Then 



\£{t,sM\<^y k>t>S. (1.14) 



Proof. It remains to check the conditions for Theorem lA.il It holds 



^ /*' ^^ ^ ^'(^)(i + ^g) _i<M. 



N 



1 + sJ, A2(r) X^{t^){l + s) -X{t^)il + .sm 
using the definition of the zone together with 



<1 



following from the monotonicity of pYrr for large t. D 

Further results for higher order derivatives. In order to perform a perfect diagonal- 
ization in the hyperbolic zone it is essential to find symbol estimates for £(t^,0,^) for 

Lemma 1.5. Assume that (Al) and (A3) holds. Then for \^\ < N and all a G N" the 
symbol like estimate 

iiD|f(t5,o,e)ii<c„^ie|-" 

is valid. 

Proof It holds BtS = AS with 

^\ . 1 



Mt,0 = (i+tu^2 it\ , ll^(*,OII 



< 



Thus for |a| = 1 we get 

BtB^S = B1{A£) = (D1A)£ + A(Df f ) 
or using Duhamel's formula together with the initial value D£<S(0, 0,^) =0 

Df£:(t,0,O= / £{t,T,C){B^A{T,0)£(T,0,OdT. 
Jo 

Now the known estimates for £{t, s,^), (|1.14|) . imply together with ||D?A(t,^)|| < 1 the 
desired statement ||D^£:(t,0,OII < t < \^\~' for (t,0 G Zd^ss{N) . 



1 Representation of Solutions 



For |a| = £ > 1 we use Leibniz rule to get similar representations containing all deriva- 
tives of order less than \a\ under the integral and use induction over i. 

The estimates 

for |ai| + |a2| < ^) formula (|1.14|) and the first statement we conclude 

l|D?£(«)ll</'^|{|-'dr<3^ir"-l 

with Proposition 11.31 1. By application of the equation itself we get estimates for time 
derivatives. This means 



^t^!nt,m<^\ 


f 1 V|.|-H 


Kl + t) 1^1 


iDfid^ier'-'"', 


\C\<N 



and together with 

(1.15) 

the statement follows. D 

This result can be reformulated as the symbol estimate 

X\t^)£it^,0,C)GS'^, 

where 

5^ = { m G C°^{W \ {0} I Va : |D"m(e)| < CalCl''"'"' } 

denotes the homogeneous symbol class of order k. Thus as consequence of Marcinkiewicz 
multiplier theorem the Fourier multiplier with symbol A^(t^)i5(tg, 0,^) maps L^ into L^ 
for all p £ (1, oo). 

Remark. For all proven results on £{t, s, (,) within the dissipative zone we do not need 
any estimates for the derivatives of the coefficient function b{t). The essential point is 
always the positivity of the coefficient and therefore the monotonicity of the auxiliary 
function A(t). 

1.2 The hyperbolic zone 

Diagonalization. We use the special symbol classes introduced by Definition ll.il Re- 
mark that it holds \S,\4'hyp,N £ •S'tvIIjO} and by Assumption (A2) also b(t)(j)hyp^N £ 
S'tvJO, 1}. For the further calculations we omit the characteristic function (j)hyp,N- 

Thus we consider 

U = {\^\u,Dtuf 
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with 



D.f/ = ^(U)f/=(j^| i%]U. 



We apply two transformations to this system. In a first step we diagonahze the homo- 
geneous principle part. After that we perform further diagonalization steps to make the 
remainder belong to a sufficiently nice symbol class. 

Step 1. We denote by M the matrix 

M=(J -^) (1.16) 

consisting of eigenvectors of the homogeneous principle part of A{t,^) with inverse 

M-'=^_\ ;). (i.iT) 

Then for U^^^ = M~^U we get the system 

DtC/(°) = {V{0 + Rit)) C/(°) (1.18) 

with 

-<«K" -i.i)' ^<"^^(! !)■ <-'' 

Step k + 1. We construct recursively the diagonalizer Nk{t,(,) of order k. Let 

k k 

Nk{t,i) = Y,N'^^\t,i), Fk{t,i) = Y,F^^\t,0 
i=o i=o 

where A^^o) = I, S(o) = R{t) and F(o) = diag B^o) = Fo(t). 

The construction goes along the following scheme. Remark that Fq is a multiple of /. 
Then we set 

F(^) =diagS(^\ 

^O'+i) = (Di - P - i?)iV,-+i - A^,+i(Di - P - F,). 

Now we prove by induction that A^'^-'^ G S]\f{—j,j} and i?^-'-' G Sn{—J,J + 1}- For j = 
we know 

f(o) g 5^{0, 1}, Ar(i) G 5;v{-l, 1}, B^'^ G 5^{-l, 2}, 
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the last one follows from the representation B^^^ = D^A^^^) — {R — F^^^)N^^\ 

For j > 1 we apply the principle of induction. Assume we know B^^' £ Sn{—J,J + !}• 
Then by definition of Ar(i+i) ^^ ^iave from |^|-^ G S'Ar{-l,0} that iV^^+i) G Sn{-J - 
l,j + 1} and F(j) G Sn{-J,J + !}• Moreover, 

^(i+i) = (D^ _ p _ ^)(^ ^H) _ (^ iv('^))(Dt - 2? - ^ F^-'')) 

1^=0 1^=0 u=Q 

+ iv(J'+i) V F^'^) - (V iv^'^) )f(j') . 

Now i?"^ + [N^^^^',!)] — F"' = for all j. The sum of the remaining terms belongs to 
the symbol class Sn{-J - 1, j + 2}. Hence B^+i) G Sn{-3 - 1,J + 2}. 

Now the definition of B^ > implies the operator identity 

{V)t-V{i)-R{t))Nk{t,i) = Nk{t,i){l}t-V{i)-Fk-i{t,0) modSjv{-fc,fc + l}. 

(1.20) 

Thus we have constructed the desired diagonalizer if we can show that the matrix Ni;{t, S,) 
is invertible on Zhyp{N) with uniformly bounded inverse. But this follows from Nk — / G 
StvJ— 1, 1} by the choice of a sufficiently large zone constant N . Indeed, we have 

IliVfc - /|| < C^6(t) < ^'^ITTt) " ° as TV ^ oo. 

Thus with the notation Rk{t,^) = —Nj^ {t,^)B^ '[t,^] we have proven the following 
lemma. 

Lemma 1.6. Assume (Al) and (A2). 

For each A; G N there exists a zone constant N and matrix valued symbols 

• Nk{t,C) G 57v{0,0} invertible for all (t,^) and with N~'^{t,^) G 5Af{0,0} 

• Ffc_i(t,0 G Sn{0, 1} diagonal with Ffc_i(t,^) - ^/ G 5Ar{-l,2} 

• Rk{t,O^SN{-k,k + l}, 

such that the (operator) identity 

(Dt - V{0 - R{t))Nk{t,0 = Nk{t,0{Bt - V{0 - Ffc_i(i,0 - Rk{t,0) (1-21) 
holds for all (t,^) G Zfiyp{N). 



10 
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Restricted smoothness assumptions. It is not necessary to know this lemma for all 

A: G N, thus we can replace Assumption (A2) by 



{A2)e for all A; = 0, 1, . . . , £ it holds 



k 



d, , ^ 



<Ck 



1 + t 



i+fc 



i.e. the corresponding one with finite smoothness. This gives b{t)4>hyp,N £ Sj^°°{0, 1} 
and the above lemma remains true under Assumption (A2)£ for all k < i and with 
matrices Nk £ S^-'^+^'^jO.O}, Fk-i G 5^-'=+^'°°{0, 1} and Rk G S^-''''^{-k,k + 1}. 

Remarks on perfect diagonalization. Lemma [0)1 can be understood as perfect diago- 
nalization of the original system. If we define F{t,^) as asymptotic sum of the F^^'{t,^), 

oo 

F(t,e)~j;F('=)(t,0, (1.22) 

fc=0 

this means we require F{t, ^) — Fk{t, ^) G SN{—k — 1, A; + 2} for all /c G N, and similarly 

oo 

iv(t,e)~J]iv('=)(t,e), (1.23) 

fc=0 

which can be chosen to be invertible, equation (|1.2fl|) implies 

{Bt - V{0 - R{t))N{t,0 - N{t,0{Bt - V{0 - F{t,0) G [] SN{-k,k + 1}. (1.24) 

fceN 

Thus if we define the residual symbol classes 

n{m}:= f] SN{mi,m2} (1.25) 

mi+m2=m 

we can find Poo(t,^) G ^{1} such that 

{Dt - V{0 - R{t))N{t,0 = N{t,0{Dt - V{0 - F{t,0 - Poc{t,0)- (1-26) 

The classes TC{m} are invariant under multiplication by exp(ibzt|,^|). This explains why 
we perform more than one diagonalization step. Multiplication by e '^' is not a well 
defined operation on the symbol classes SN{mi,m2}, it destroys the symbol estimates 
according to the following proposition. 

Proposition 1.7. 

1. e±**l«l5^'^2{mi,m2}-^5^'^2{mi+£,m2-^} with I = £i + £2 
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2. e±**l«lW{m} ^ n{m\ 
Proof. It suffices to prove tlie first statement. It holds for a G Sj\[{m'i,m2} 

ki+k2=k ai+a2=a 

/ 1 \ m2+fc2-«i 

< \^ \^ C' |t|mi-a2+fci / _ 

— Z^ Z-/ ^'fci,fc2,ai,a2l''l \l + t 

fci+fc2=fe ai+a2=a 

1 \ m2-i-\-k 



for k < ii, a < £2 using Leibniz rule and the definition of the hyperbolic zone. D 

Fundamental solution of the diagonalized system. After performing several diagonal- 
ization steps we want to construct the fundamental solution of the diagonalized system 

{Bt-V{O-Fk-i{t,O-Rk{t,O)Sk{t,s,O=0, ffc(s,s,e)=/GC2>^2 (;l_27) 
and obtain structural properties of it. The construction goes along the following steps: 

• the fundamental solution £o{tj s, ^) to Dj — P(C), 

• influence of the main term F^^'(t, ^) of Ffe_i(t, ^), 

• influence of Fk{t,0 - ^^"H^^O and Rkit,0- 

The fundamental solution £Q{t,s,S,) describes a 'phase function' of a Fourier integral 
operator, i.e. the oscillatory behaviour of the solution multiplier. The main term 
-^ (^lO describes the energy decay. Together with the other terms it constitutes a 
Fourier multiplier which behaves as symbol with restricted smoothness. The number k 
of diagonalization steps stands in direct connection to the smoothness properties of this 
symbol. 



Step 1. Let 



i(t-s)|C| 



£oit,s,0=exp{iit-s)ViC)} = ( ^-^it-sM\ ) (1-28) 



such that for £Q{t,s,^) = -^vSoitjS,^) the equation 

Bt£o{t,s,0 = {V{()+F^^\t,0)£o{t,s,0 (1.29) 

is satisfied. Thus Sq describes the influence of the main diagonal term. 
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Step 2. By the aid of £o{t,s, $,) we define 

nk(,t,s,o = £o(.s,t,0{Fk~i(.t,0 + Rk{t,0 - F^''Ht,0)So{t,s,o, 

= Fk-i{t,0 +£o{s,t,^)Rk{t,0£o{t,s,0 - F'^'>\t,0, (1.30) 

such that using the solution to 

BtQk{t,s,0=nk{t,s,OQk{t,s,0, Qfc(s,s,0 = /GC2>^2 (1.31) 

the matrix <Sfc(t, s, ^) can be represented as 

£k{t,s,0 =£o{t,s,OQk(.t,s,0 = j^£o{t,s,OQk{t,s,0- (1.32) 

The solution to ()1.31() is given by the Peano-Baker formula as 

°° ft fti i'te_i 



Qk{t,s,0=I + Y.^^ I ^k{ti,s,0 rnk{t2,s,()... r 'nk{te,s,Odte...dti. 

(1.33) 
To estimate this representation we use the following well-known inequality. It holds 



Proposition 1.8. Assume r G Lj^^(M.). Then 

r{ti) I r{t2)... r \{tk)dtk...dti 

J s J s 

for all A; G N. 



<^(^*Kr)|dr) (1.34) 



Step 3. The series representation (|1.33() for Qk(t, s, S,) can be used to deduce estimates. 
From the unitarity of £o{t, s, ^) it follows that 

\\nk{t,s,o\\ = \\Rkit,0\\< ^ 



and thus using 

°° dr _ 1 _ 1 

,^ {l + rm\- (l+t^m-N 
from the representation of Qk by Peano-Baker formula it follows 

\\Qkit,s,m<^- 

In a second step we want to estimate ^-derivatives of Qk(t,s,£,). Proposition II .71 yields 
for i?fc(t,^) G Sj^ '^{—k,k + 1} under Assumption (A2)^ with 

k-l<£-k 

that TZk{t,s,^) £ S^ ' ~ { — 1,2} uniform in s (and derivations w.r.to s behave like 
multiplications by |^|). 
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Proposition 1.9. Assume a € 5^ {—1, 2}. Then 



'N 
°° ft ft 



"^ f'T f'Ll ft ' 1 

b{t,s,^) = l + Y, a{ti,0 a{t2,0---l a{tj,Odtj...dti 

•1 J s J s J s 

defines a symbol from Sj^ {0, 0} uniform in s > t^. 

Proof. We use Proposition II . 81 to estimate this series. This yields in a first step (without 
taking derivatives) 

and taking a derivatives with respect to ^ yield in each summand further factors |^|~'°' 
according to Leibniz rule. D 

An almost immediate consequence of this proposition is the following structural repre- 
sentation of the fundamental solution. 

Theorem 1.10. Assume (Al) and (A2)2k-i, k > 1. Then the fundamental solution 
^kitiS,^) of the diagonalized system H1.27|) can be represented as 

X(s) 

£k{t,s,0 = j7w£o{t,s,OQk{t,s,0 t,S>t^ 

with a symbol Qkit,s,^) of restricted smoothness subject to the symbol estimates 



BlBlQk{t,s,0 <C,,,|Crl"l t,s>t^ 

for all multi-indices \a\ < k — 1 and all ^ G Nq. 

The first derivative with respect to t can be estimated by the equation for Qk directly. 

Of special interest is £kit,t^,£,). The estimate of the previous lemma together with the 
properties of the derivatives of t^ imply 

Corollary 1.11. Assume (Al) and (A2)2k-i, k>l. Then 

Qk{t,t^,O^S'/~'{0,0} 

(for t > t/: and \C\ < Nf 

^We could also make tj a smooth function in ^ to avoid such problems. But then we can not think 
about fundamental solutions £{t,s,^) with variable starting time a. 
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The matrix Qfcl^;*,^) converges for i ^> co to a well-defined limit. This limit will be 
used in Section IHl to conclude the sharpness of our results. 

Theorem 1.12. Assume (Al) and (A2)2k~i, k>l. The limit 

Qk{oo, s, = lim Qk{t, s, C) 

exists uniform, in ^ for \^\ > ^g. Furthermore 

||DfQfc(oo,tg,OII<C,|e|-l"l 
for all multi-indices \a\ < k — 1. 

Proof. We fix the starting value s and consider only |^| > ^g (i.e. s > t^). Taking the 
difference Qk{t,s,^) — Qk(t,s,^) in the series representation (|1.3.Sj) yields 

oo „t 



Qkit,s,^)-Qkii,s,0 = Yl I ^kiti,s,0 r nk{t2,s,i)... r 'nkite,s,C)dti...dti 

■_i Ji Js J s 

such that with Proposition II. SI 

\\Qkit,s,0 - Qk{i,s,mL^m>^.} < J, ||i?(ti,e)||expM~||/?(T,OI|drldii 



fOO 

0, t,i ^ oo. 



Similarly one obtains for \a\ < k — 1 

||D^Qfc(t,s,0-D^Qfc(t,s,OII<l^r"^ ^IT^^O t,t^oo 
uniform in |,^| > ^^ Now the second statement follows from the estimates of t^, formula 

We have even proved more. The limit exists in the symbol class S^ of restricted smooth- 
ness k — 1. 

Proposition 11.81 mav also be used to estimate the formal representation of Qk{oo,s,S,) 
as symbol in (s,^). 

Corollary 1.13. The series representation 



/'OO /•CI ri^£~i 

Qk{oo,s,0 = I + ^i^ 7^fe(tl,s,0/ nkit2,s,0--- nkiti,s,Odti... 

j^i-^S 'J S 'J s 



dti 



^0,fe-l 



gives an asymptotic expansion of Qj,{oo,s,S,) in Sj^ {0,0}, i.e. the jth term belongs 
to S^/~^{-j,j}. 
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2 Estimates 



Step 4. The transpose of the inverse of Qk satisfies the related equation 

DtQfit,s,C)+nl{t,s,OQf{t,s,0=0, Qf{s,s,0=l€^C'''\ (1.35) 

The matrix TZ^{t,s,^) satisfies the same estimates like TZkitjS,^) and therefore the 
reasoning of the previous step holds in the same way for Q^ (t,s,^). Especially the 
matrix Qk{t,s,^) is invertible in the hyperbolic zone and Q^ {oo,s,^) exists. 

Corollary 1.14. Assume (Al) and (A2)2k-i, k'>l. Then the limit 

exists uniform, in ^ for \^\ > S,s- 

Transforming back to the original problem. After constructing the fundamental solu- 
tion £k{tj s, ^) we transform back to the original problem and get in the hyperbolic zone 
the representation 

£{t,s,0 = MNk{t,0£k{t,s,ON^\s,OM-^ (1.36) 

with uniformly bounded coefficient matrices Nk,N^ G S'tvJO, 0}. We combine this 
representation with the representation obtained in the dissipative zone. This yields for 
s < t^ the representation 

£{t,s,0 = -^^MNk{t,C)£o{t,t^,OQit,t^,mk'ik^OM~'X{t^)£{t^,s,C). (1.37) 

2 Estimates 

The so far obtained representations of solutions allow us to conclude estimates for the 
asymptotic behaviour. This section is devoted to the study of estimates which are 
directly related to our micro-energy ()1.6|) . i.e. estimates for the fundamental solution 
£{t,s,D) or to the closely related energy operator E(i,D). 

2.1 L'^-L'^ estimates 

Theorem 11.101 together with the estimate of Lemma ll .41 implies by Plancherel's theorem 
the following operator norm estimate. 

Theorem 2.1. Assume (Al), (A2)i and (A3). Then the L?-L? estimate 

ll«(M.D)|b..<^ 
holds. 
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2 Estimates 



Using the definition of the micro-energy ()1.6|) we can reformulate this estimate in terms 
of the energy operator E(i, D). For convenience we recah the relation between the 
multiplier £{t,s,S,) and E(t,^). They are a direct consequence of the definition of our 
micro-energy (|1.6|) . 

Proposition 2.2. 1. It holds E(f , ^ = £{t, s, OH^, for s>t^. 
2. 

3. The multiplier |C|//i(t, ^) induces a uniformly bounded family of operators on L^, 
p G (1, co) converging strongly to the identity for t — > co. 

Corollary 2.3. Assumptions (Al) (A2)i, (A3) imply 

< 1 



l|E(t,D)||2^2^ 



m 



We conclude this section with examples. 
Example 2.1. Let 

Then Assumptions (Al) (A2) and (A3) are satisfied and the above Corollary gives again 
the known estimate ||E(t,D)||^^2 ^ (1 + *)~^ from |Wirn4j . 
Example 2.2. Let fi > and n > 1. Then we consider 

Ht) = r^-i • 

(1 + t) log(e + *)••• log["] (eN + t) 
Again the assumptions are satisfied and we obtain 

A(t) = (logW(eM+i))2 

and the energy decay rate may become arbitrary small. This example is taken from the 
paper of K. Mochizuki and H. Nakazawa, |MN96j . 

Example 2.3. To give at least one example with an oscillating coefficient we consider 

^ 2-Fcos(alog(e-Ft)) 
^'~ 4(e + t) 

with large real a. Then it holds 

/■ 1 1 

/ b{t)dt = - log(e + t) + -— sin(a log(e + t)) 
J 2 4a 



and thus 



X(t) = Cq \/ e -|- t exp < - — sin(alog(e -|- t)) 
[ 4a 



such that the energy decay is given by ||E(t,D)||2^2 ~ (e + 1) 2^ which is independent 
on the choice of a. 
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2 Estimates 



2.2 LP-L"? estimates 

This section is devoted to L^-L'^ estimates. The basic estimate is given in the following 
theorem, it restates the known result for free waves in the language of our operators. 

Theorem 2.4. It holds 

||fo(i,D)||p,,^g<Cp,,(l + t)"'^(p~i) 

for dual indices p and q, p G (1,2] and with regularity r = n ( 

By the aid of this estimate we deduce from our representation a corresponding estimate 
for the dissipative equation. 

Theorem 2.5. Assume (Al), (A2) and (A3). Then the operator £{t,s,D) satisfies for 
dual indices p and q, p & (1, '^], pq = P + q the norm estimate 

with regularity r = n (- — - 

Proof. We decompose the proof into two parts. First we consider £{t, 0, 1))(l)diss,N{t, D)- 
Using the estimate \\£{t,0,^)cl)diss,N{'tiO\\ ^ wiy together with the definition of the 
zone we get 

which is a stronger decay rate than the one given in the theorem. 

In a second step we consider the hyperbolic part. For small frequencies we use the 
representation 

£{t,^,C)4>hyp,N{t,i) = 

-i- MNk{t, e) £o{t - t^, Qk{t, %, N~\t^,OM-^X{t^)£it^,0, <PhypAi^ 

A(t) > ^ 'V ^ 'V ^ ^v ^ ' 

q—*q 'p,'i — >q p,'i — >p,r p,r—>p,r 

together with the mapping properties of the multipliers marked with a brace. It is 
essential that k — 1 > [^] . The operator £o{t, ^) brings the hyperbolic decay rate, while 

For large frequencies the representation simplifies to 

£{t,0,0<Phyp,N{t,0 = tL MNk{t,0£o{t,^) Qk{t,0,0 N~\0,^)M-^ cl)hyp,Nit^O, 

^[^) '• V ■'^ V '^ V ■'^ V ' 

g— >g p,r-^q p,r^p,r p,r^p,r 

the argumentation remains the same. D 
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2 Estimates 



Minimal regularity for the W-L'^ estimate. With the notation 

^ I n + 1, n odd, 

we can prove the above U'-L'^ decay estimate under the weaker Assumption {A2)i^ on 
the coefficient function. If we use this regularity of the coefficient and perform k = [^] 
diagonahzation steps we obtain Nk{t^^) G 5'j^°°{0,0} and Qk{t-,s,S) is uniformly in 
t > s > t^ a symbol of smoothness [^] . Thus, Nk(t>hypj N^ (j)hyp and Qk4>hyp define 
operators L^ -^ L^ for all p G (1, oo) with uniformly bounded operator norm in t, s. 

2.3 Estimates of the solution itself 

The definition of the micro-energy allows to extract also estimates for the solution from 
the structure oi £{t,s,$,). Similar to the reasoning in the previous theorems one obtains 
for the solution operator 

normalised in such a way that S(t,D) : L^ — > L^, the following L^-L'^ estimate. 

Theorem 2.6. 1. Assume (Al), (A2)i and (A3). Then the solution u(t,x) satisfies 

the L^-L^ estimate 

l + t 



\S{t,B)\\2^2< 



x^ty 



2. Assume (Al), (A2)i^ and (A3). Then the solution u{t,x) satisfies the L'^-L'^ 
estimate 

faTY(l + i)^""(^"^), p>p\ 

[^(1+t) 2 y-^ ~^), p<p\ 

for dual indices p G (1, 2], pq = p + q and with regularity r = ""- ( p ~ ^ ) • The 

critical value p* = 2n+2 '^■f ^^(0 ~^ ^'^^ determined by the balance of the two 
estimates in the remaining cases. 

Proof. We start with the L^-L^ estimate. 

In the hyperbolic zone and for |^| > c we estimated |^|m by A~"^(t) and we can just divide 
by 1^1 to get an estimate for u by 

1 ^ l + t 



A(t)-A2(t)- 

For small |^| we have to take into account that dividing by h{t, ^) brings a further factor 
(1 + t). Thus we obtain inside the dissipative zone (1 + t)/}?(t) and in the hyperbolic 
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3 Sharpness 



zone for small frequencies 



1 + tf < 1 + ^ 



using the monotonicity of t/\{t) for large t following from (A3). 

To conclude the L^-L'^ estimate we follow the same proof like in the case for the energy 
with the exception that the rate inside the dissipative zone is different and dominating 
at least if p and q are near to two. D 

This result coincides for the case b{t) = j^, fi < 1 with the estimate from |Wirfl4j . 
We think the following observation is important to understand the essential difference 
between the solution estimate and the estimate for the energy. 

For p near to 2 the estimate for the solution operator S(t, D) comes from 
properties of the dissipative zone, while estimates for the energy E(i, D) are 
determined by the hyperbolic zone. 

Remark 2.4. Assumptions on the data, which are effective near the frequency ^ = 
may be used to improve these estimates. One possibility to achieve this is the use of 
weighted initial data. If we assume that (x) ui,{x) U2 S L^ and the space dimension 
n > 3, then the estimate inside the dissipative zone in this theorem might be improved 
by an application of Hardy inequality and one obtains 

3 Sharpness 

Finally we want to prove the sharpness of the above given energy decay estimates. Our 
constructive approach enables us to formulate the question of sharpness as a modified 
scattering theory. The basic idea is as follows 

• we relate the energy operator E{t,D) to the corresponding operator Eo(t,D) for 
free waves multiplied by the decay rate 

• this relation defines a M0ller wave operator W+ defining appropriate data to the 
free wave equation with the same asymptotic properties (up to the factor A(i)), 

• furthermore, we need to know the mapping properties of the M0ller wave operator, 

• and the convergence defining the wave operator has to be understood. 

A first observation follows immediately from Liouville Theorem and gives an expression 
for the determinant of E(t, ^). 
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3 Sharpness 



Lemma 3.1. It holds detE(t,^) = j^iC], where [£,] = \C\/ (0- 

After these introductory remarks we can state the following theorem. It holds 
Theorem 3.2. Assume (Al), (A2) and (A3). Then the limit 

VF+(D) = s-limA(i)(Eo(t,D))-iE(t,D) 

t— >oo 

exists as strong limit in Lp' -^ L? and defines the modified M0ller wave operator W+. It 
satisfies 

w+io = {Eoit^,or'Qkioo,t^,ON^Hk,OM~'x{t^Mk,c) 

for all k > 1. 

Remark that t^ depends on the zone constant and this constant is chosen after diago- 
nalizing k steps. Thus 

Qk{^,t^,ONk\k,0 

is independent on 1 < A; < ^ for sufficiently large zone constant A^ depending on i. 
Proof. The proof consists of three steps. 

Step 1. With the notation 

Vc = {U e L^\ dist(0, supp U) > c} 

we can construct the dense subspace M = IJc>o^c- Now Theorem 11.121 together with 
the representation E(t,^) = £{t,t^,^)E,{t^,^) implies the existence of the limit 

lim A(t)Eo(t,D)~^E(t,D) 

t— >oo 

as limit in the operator norm in Vc ^ Vc for all c > 0. Thus the limit exists pointwise 
on M. 

Step 2. The energy estimate, Theorem 12.31 implies that A(t)Eo(i,D)^^E(t,D) is uni- 
formly bounded in I? ^f I?' . Thus the Theorem of Banach-Steinhaus implies the exis- 
tence of the strong limit and defines VF+. 

Step 3. The previously defined operator VF+ is given on each subspace Vc as Fourier 
multiplier with symbol 

w+{i) = {Eoit^,or'Qk{^,k,ON^Hk^OM-'\{t^)E{t^,o, 

which is independent of c. Thus the representation holds on M and using the bounded- 
ness of W4. on the whole space. D 
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A Appendix: A remark on Volterra integral equations 



A(t)£(t,0,D) 



So{t,D) 




t = 



Figure 3.1: Modified scattering theory. 



Remark that for above given theorem we need (A2)2fc-i to represent Ty+(^) in terms of 
Qkit, s, (,). Thus Assumption (A2)i is sufficient to get the existence of the wave operator. 

Corollary 3.3. It holds 

detW+{0 = l, 

and therefore W+(D) defines an isomorphism on L^ . 

Interpretation of the result. Theorem 13.21 mav be used to construct for each datum 
{{D) ui,U2) G L^ to Cauchy problem (|0.1|) a corresponding datum ((D) ■ui,-U2)"^ = 
W+(D){{D) ui,U2)'^ to the free wave equation Du = 0, such that the solutions are 
asymptotically equivalent up to the decay factor A~^(t), i.e. it holds 



|Eo(i, D)((D) ui,ii2f - A(t)E(t, D)((D) m, U2 



iTI 







as t 



c«. 



(3.1) 



This is a direct consequence of the unitarity of Eo(t,^). It implies that the above given 
L^—L^ estimates are indeed sharp and describe for all nonzero initial data the exact 
decay rate. 



A Appendix: A remark on Volterra integral equations 



We are interested in solutions to the Volterra equation 



f{t,p)=f{0,p)+ / k{t,T,p)f{T,p)dT 

Jo 



(A.l) 



with kernel k = k{t,T,p) and given initial value f{0,p) depending on some parameter 
pe PC M". 
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References 

Theorem A.l. Assume /(O, •) G L°°{P), k e L°°{R\ x P) and 

I |A;(t,r,p)|dr G L°°(M+ x P). 

Then there exists a (unique) solution f{t,p) of ()A.1|) in L°°(M+ x P). 

For the proof of this statement we follow |GLS9nj . The condition on the kernel implies 
that the operator 

K : f{t,p)^ [ k{t,T,p)f{T,p)dT 

Jo 
is bounded. This implies that, if the norm of this operator is sufficiently small, the 
solution of the integral equation is given by the contraction mapping principle. Now 
the idea of |(lLS9n| Chapter 9.3, Theorem 3.13] applies, we can decompose the time 
interval into a finite number of smaller sub-intervals where the restricted operators are 
contractions. This follows from 



esssup [ ' \k{t,T,p)\dT < (T2 - Ti) ||A;||oo ^ 0, as T2 



Furthermore, the Volterra structure of the equation can be used to build up the resolvent 
from the resolvents of the restricted operators like in |GLS90j . 
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